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Abstract: In this work, we offer a systematic framework to describe the gluodynam- 
ics and chiral dynamics. We firstly construct a quenched dynamical holographic QCD 
(hQCD) model in the graviton-dilaton framework for the pure gluon system, then develop 
a dynamical hQCD model for the two flavor system in the graviton-dilaton-scalar frame- 
work by adding light flavors on the gluodynamical background. The dimension-2 dilaton 
background field ^ takes the same form as that in the soft-wall model and its square 
can be identified as the dual field to the gauge invariant operator of the dimension-4 gluon 
condensate. The gluon condensate in the QCD vacuum breaks conformality and induces 
a deformed warp factor. Without introducing extra parameters but just self-consistently 
solving the deformed metric induced by the dilaton background field, we find that the 
scalar glueball spectra in the quenched dynamical model is in very well agreement with 
lattice data. For two flavor system in the graviton-dilaton-scalar framework, the deformed 
metric is self-consistently solved by considering both the chiral condensate and gluon con- 
densate in the vacuum, which are responsible for the chiral symmetry breaking and linear 
confinement, respectively. It is found that the mixing between the chiral condensate and 
gluon condensate is important to produce the correct light fiavor meson spectra. The pion 
form factor and the vector couplings are also investigated in the dynamical hQCD model. 
Besides, we give the criteria for the existence of linear quark potential from the metric 
structure, and show a negative quadratic dilaton background field is not favored in the 
graviton-dilaton framework. 
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1 Introduction 

The local quantum field theoretical description has made great success since it was firstly 
developed in the quantization of electrodynamics and further been developed and applied 
to the descriptions of elementary particles. Nowadays, quantum chromodynamics (QCD) 
is accepted as the fundamental theory of the strong interaction. In the ultraviolet (UV) or 
weak coupling regime of QCD, the perturbative calculations agree well with experiment. 
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However, in the infrared (IR) regime, the description of QCD vacuum as weh as hadron 
properties and processes in terms of quark and gluon stih remains as outstanding challenge 
in the formulation of QCD as a local quantum field theory. In order to derive the low- 
energy hadron physics and understand the deep-infrared sector of QCD from first principle, 
various non-perturbative methods have been employed, in particular lattice QCD [1-4], 
Dyson- Schwinger equations (DSEs)[5, 6], and functional renormalization group equations 
(FRGs)[7-9]. 

In recent decades, an entirely new method based on the anti-de Sitter /conformal field 
theory (AdS/CFT) correspondence and the conjecture of the gravity/gauge duality [10-12] 
provides a revolutionary method to tackle the problem of strongly coupled gauge theories. 
Though the original discovery of holographic duality requires supersymmetry and con- 
formality, the holographic duality has been widely used in investigating hadron physics, 
strongly coupled quark gluon plasma and condensed matter. It is widely believed that the 
duality between the quantum field theory and quantum gravity is an unproven but true 
fact. In general, holography relates quantum field theory (QFT) in d-dimensions to quan- 
tum gravity in (d + l)-dimensions, with the gravitational description becoming classical 
when the QFT is strongly-coupled. The extra dimension can be interpreted as an energy 
scale or renormalization group (RG) flow in the QFT [13-20]. 

Many efforts have been invested in applying holography duality to describe the real 
QCD world, e.g. for mesons [21-23], baryons [24] and glueballs [25], see Refs. [26, 27] for 
reviews. It is well-known that the QCD vacuum is characterized by spontaneous chiral 
symmetry breaking and color charge confinement. The chiral symmetry breaking can be 
read from the mass difference between the chiral partners of the hadron spectra, and 
the spontaneous chiral symmetry breaking is well understood by the dimension-3 quark 
condensate {qq) [28] in the vacuum. The understanding of confinement remains a challenge 
[29]. From the hadron spectra, confinement can be read from the Regge trajectories of 
hadrons [30, 31], which suggests that the color charge can form the string-like structure 
inside hadrons. Confinement can also manifest itself by the linear potential between two 
quarks at large distances [32]. 

A successful holographic QCD model should describe chiral symmetry breaking, and 
at the same time should describe both the Regge trajectories of hadron spectra and linear 
quark potential, two aspects in the manifestation of color confinement. Nonetheless, these 
important nonperturbative features haven't been successfully accommodated in a unique 
hQCD model. 

The current achievements of AdS/QCD models for hadron spectra are the hard- wall 
AdS/QCD model [21] and the soft-waU AdS/QCD or KKSS model [22]. In the hard- wall 
model [22], the chiral symmetry breaking can be realized by chiral condensation in the 
vacuum, however, the resulting mass spectra for the excited mesons behave as ~ n^, 
which is different from the linear Regge behavior ~ n. In order to generate the linear 
Regge behavior, the authors of Ref.[22] introduced a quadratic dilaton background, one can 
obtain the desired mass spectra for the excited vector mesons, while the chiral symmetry 
breaking phenomenon cannot be consistently realized [33]. 

Interesting progress was made in Refs. [34, 35], where a quartic interaction term in the 
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bulk scalar potential was introduced to incorporate linear trajectories and chiral symmetry 
breaking. However, such a term was shown [35] to result in a negative mass square for the 
lowest lying scalar meson state, which might cause an instability of the background. In 
Ref.[36], a deformed warp factor is introduced, which can cure the instability and maintain 
the linear behavior of the spectra. 

With AdSs metric in the soft- wall model and its extended versions [34-41] (except [36]), 
only Coulomb potential between the two quarks can be produced [42]. On the other hand, 
the linear quark potential can be realized in the Andreev-Zakharov model [43], where a 
positive quadratic correction in the deformed warp factor of AdSs geometry was introduced. 
The linear heavy quark potential can also be obtained by introducing other deformed warp 
factors as in Refs. [44, 45]. The positive quadratic correction in the deformed warp factor 
in some sense behaves as a negative dilaton background in the 5D action, which motivates 
the proposal of the negative dilaton soft-wall model [46, 47]. More discussions on the sign 
of the dilaton correction can be found in [48, 49]. 

It is noticed that the quadratic correction, whether appears in the 5D action or in the 
deformed warp factor, indeed plays an important role to realize the linear confinement, 
though only partly. Since both the Regge trajectories of hadron spectra and linear quark 
potential are two aspects in the manifestation of color confinement, they should share the 
same dynamical origin and should be realized in the same holographic QCD model. 

In the soft- wall model and its improved versions, the dilaton field or the deformed 
warp factor are introduced by hand. In Ref.[50], we have successfully described the chiral 
symmetry breaking, the Regge trajectories of hadron spectra and linear quark potential in 
the graviton-dilaton-scalar coupling framework, in which the metric, the field(s) and the 
potential(s) of the field(s) are self-consistently determined by field equations, and one can 
self-consistently solve out the other two with one input. There are at least three different 
ways to deal with the system in the literature: 1) Input the form of the field(s) to solve 
the metric structure and the potential(s) of the field(s) [50, 51]; 2) Input the potential(s) 
of the field(s) to solve the metric and the field(s) [52, 53]; 3) Input the form of the metric 
structure to solve the field(s) and the potential(s) of the field(s) [54]. 

This work is an extension of Ref . [50] , and we intend to establish the relation between 
the QCD dynamics including at IR and its induced geometry. The paper is organized as 
follows: In Sec. 2, we establish a quenched dynamical hQCD model in the graviton-dilaton 
framework to describe the pure gluon system, and by selfconsistently solve the deformed 
warp factor induced by the dilaton field, we get the scalar glueball spectra; We introduce 
the meson spectra in the KKSS model and its improved versions in Sec. 3; We then develop 
the graviton-dilaton-scalar coupling framework for two flavor system and investigate the 
hadron spectra in Sec. 4, and we also investigate pion form factors and vector couplings in 
Sec. 5. We give discussion and summary in Sec. 6. 
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2 Pure gluon system: quenched dynamical holographic model 

At the classical level, QCD is a scale invariant theory, which is broken by quantum fluctu- 
ations. The pure gluon part of QCD Lagrangian in 4-dimension is described by 

= -Ig%{x)G^''''^{x), (2.1) 

with 

G;;,(x) = d.A^ix) - d,A;{x) + gr'"^A'^{x)Al{x). (2.2) 

Where Afj^{x) is the gluon field with a = 1, • • • ,8 the color indices. 

In the vacuum, the scale invariance is explicitly broken, and the relevant degrees of 
freedom of QCD at infrared are still poorly understood. Varies of vacuum condensates 
provide important information to understand the non-perturbative dynamics of QCD. For 
example, the gauge invariant dimension-4 gluon condensate {g^G"^) has been widely investi- 
gated in both QCD sum rules and lattice calculations [55-57], and the non- vanishing value 
of the condensate does not signal the breaking of any symmetry directly, but rather the 
non-perturbative dynamics of strongly interacting gluon fields. In last decade, there have 
been growing interests in dimension-2 gluon condensates {g^A^) in SU(A'^c) gauge theory 
and its possible relation to confinement [58-74]. 

On the other hand, the effective Lagrangian for pure gluon system can also be con- 
structed in terms of lightest glueball [75-77] or one scalar particle - dilaton [78-83]. The 
dilaton field is an hypothetical scalar particle predicted by string theory and Kaluza-Klein 
type theories, and its expectation value probes the strength of the gauge coupling. In 
Ref.[80], an effective coupling of a massive dilaton to the 4-dimensional gauge fields pro- 
vides an interesting mechanism which accommodates both the Coulomb and confining 
potentials between heavy quarks. In the following, we construct a 5-dimension dynamical 
hQCD model in the graviton-dilaton coupled system for the pure gluodynamics. 

2.1 Quenched dynamical holographic model and gluodynamics 

Csaki and Reece in Ref.[51] proposed to model the pure gluon system in the graviton- 
dilaton framework by considering the correspondence between the dilaton background field 
and the non-perturbative gluon condensate, which provides a natural IR cut-off. They 
investigated the case of dilaton field coupling with dimension-4 gluon operator Tr(G^) and 
higher dimension-6 gluon condensation Tr(G^), and found that such IR correction cannot 
generate Regge spectra of glueball. They also discussed a tachyon-dilaton-graviton system, 
where the tachyon corresponds to a dimension-2 gluon condensate {g'^A^) to realize the 
linear confinement. However, the local dimension-2 gluon operator is not gauge invariant. 

We construct a 5D effective model for pure gluon system by introducing one scalar 
dilaton field ^{z) in the bulk. It is not known how the dilaton field should couple with the 
gauge field in 4-dimension. Motivated by Refs.[77] and [80], we introduce the holography 
dictionary as <I>^(z) dual to the gauge invariant dimension-4 gluon condensation Tr(G^). 
In this case, though the dilaton field $(z) itself has dimension of 2, the action is always in 
terms of <I>^ thus there is no gauge invariant problem. 
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The 5D graviton-dilaton coupled action in the string frame is given below: 

Sg = / d^x^se^^"" {Rs + AdM^d^'<^ - VS{<^)) . (2.3) 

Where G5 is the 5D Newton constant, Qs, $ and Vq are the 5D metric, the dilaton field 
and dilaton potential in the string frame, respectively. The metric ansatz is often chosen 
to be 

aim = bl{z){dz'' + r?^,dx^dx"), 6,(z) = e^»(^). (2.4) 

In this paper, the capital letters like "M,N" would stand for all the coordinates(0,l,..,4), and 
the greek indexes would stand for the 4D coordinates(0,...,3). We would use the convention 

rp^ = mo = -1,?/*^ = = ^ij-) 

Under the frame transformation 

9mn — 9mn'^ i — ^ ^G^ l^-^J 

Eq.(2.3) becomes 

SE = / d'x^ [re - \d^^d^^ - Vii^)^ . (2.6) 

The equations of motion can be easily derived by doing functional variation with respective 
to the corresponding fields. It takes the familiar form in the Einstein frame, 

Emn + \gL + If ($) j - ^dm^dn^ = 0, (2.7) 

and 

.r^dmiVd^d^^) - d^Vii^) = 0. (2.8) 
3V5E 

The simplest dimension-2 dilaton background field has the quadratic form as 

$ = ±fiGz'^. (2.9) 

The positive quadratic dilaton background is the same as the one introduced in the KKSS 
model [22]. We will show in Section 2.2 and 2.3 that only positive quadratic correction can 
generate the linear confinement. 

Under the metric ansatz Eq.(2.4), the above Einstein equations has two independent 
equations, 

-AE + A'i-U'' = 0, (2.10) 

+ 3Ae<^' - le^^^a^lf ($) = 0. (2.11) 
o 

in the new variables of 

bE{z) = fes(2)e-i*(^) = e^^(^), Ae{z) = As{z) - h{z). (2.12) 
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In the string frame, the above two equations of motion become 

- X', - ^^'X, + Al + = 0, (2.13) 

+ (3^1 - 2«>>' - ^e2^^-t*9<i,(et*Vg($)) = 0. (2.14) 
8 

With the quadratic dialton background given in Eq.(2.9), we can solve the metric Ae 
and the dilaton potential Vq{'^) in the Einstein frame as 

Ae{z) = log(^) - log(oFi(5/4, ^)), (2.15) 

, 12oFi(l/4,^)2 16oFi(5/4,^)2$2 
Vim = + ' , (2-16) 

with 0-^1 (a; z) the hypergeometric function. It is noticed that in the Einstein frame, both 
the positive and negative quadratic correction give the same results. However, in the string 
frame, the positive quadratic dilaton background field = ^IqZ^ gives: 

At = Ae{z) + \^?az^ Vt = e-'l^^l^\B. (2-17) 
and the negative dilaton background field = —fi^z^ gives 

AJ = Ae{z) - \nlz\ y^'- = e^/3^c^V#. (2.18) 

We will show in Section 2.2 and 2.3 that positive and negative quadratic dilaton background 
will induce different results on the glueball spectra and the quark-antiquark potential. 
With the normalized variable $Ar which is defined as 

3- '-A/^ . ^a..if.p>M. 



--dM^d'''^^--dM^d''''^, (2.19) 



the dilaton potential in the Einstein frame takes the form of 



Vi{^) = -12 + 2 0-^1 (V4, g)^^^ ^ ^2.20) 



oFi(l/4^ oFi(5/4;^ )2$ 

here L the radius of AdS5 and ^Fi^a^z) the hypergeometric function. In the ultraviolet 
limit, 

Vi + \mI^\ (2.21) 
with the 5D mass for the dilaton field 

M|l2 = -4. (2.22) 

From the AdS/CFT dictionary A(A - 4) = MgL^, one can derive its dimension A = 2, 
which is consistent with our ansatz that $2 corresponds to dimension-4 gluon condensate. 
It is worth to emphasize that even though the dilaton field itself has dimension of 2, 
however, because it always appears in the action in the form of there is no gauge 
invariant problem in our framework. 

We call this 5D model for pure gluon system as the quenched dynamical soft-wall 
holographic model. The gluon content of the dimension-2 dilaton field itself is not clear to 
us yet, there are several possibilities: 
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^ /'^^-+^ 


L/ati 






ljat4 Ijato 




iVc = 3 


Ne = 3 


Nc ^ oo 


Nc = 3 iVc = 3 


1 


1475(30) (65) 


1580(11) 


1480(07) 


1730(50)(80) 1710(50)(80) 


2 


2755(70)(120) 


2750(35) 


2830(22) 


2670(180)(130) 


3 


3370(100)(150) 








4 


3990(210)(180) 









Table 1. Lattice data for 0'^'^ glueball in unit of MeV. Latl data from Ref.[91], Lat2 and Lat3 
data from Ref. [92], Lat4 [93] and Lat5 [94] arc anisotropic results. 



1. The dimension-2 dilaton might be dual to the dimension-2 gluon condensate {g^A^) 
[61, 64, 65], which has been discussed in some literatures, e.g. Refs.[35, 39, 43, 51]. 
However, it encounters the gauge invariant problem. 

2. The simplest gauge invariant dimension-2 gluon operator is the zero momentum mode 
of {g^A^), i.e. < g^A^{k = 0) >, the Bose-Einstein condensation (BEC) of the 
"pairing" of two gluons in the vacuum due to the strong interaction [58, 68]. The 
BEC of the " pairing" of two gluons spontaneously generates an effective gluon mass 
and breaks scale invariance, and in this scenario, the dimension-4 gluon condensation 
is proportional to the dimension-2 gluon condensation. Recent lattice results support 
a gluon mass at IR [84-86] which was proposed by Cornwall in 1981 [87] and recently 
developed in [88]. 

3. The dimension-2 dilaton field might also correspond to the gauge invariant but non- 
local operator related to topological defects in the QCD vacuum [64]. 

We leave the fundamental problem on the basic degree of freedom of QCD vacuum 
and the mechanism of confinement for further studies. In the following, we show that the 
results of glueball spectra and linear potential based on the simplest ansatz of dimension-2 
dilaton background field in Eq.(2.9). 

2.2 Scalar glueball spectrum in quenched dynamical holographic model 

The glueball spectrum has attracted much attention more than three decades [89]. The 
study of particles like glueballs where the gauge field plays a more important dynamical 
role than that in the standard hadrons, offers a good opportunity of understanding the 
nonperturbative aspects of QCD, e.g. see reviews [90]. In Table 1, we list the scalar 
glueball spectra from several lattice groups [91-94]. 

The glueball has been studied in the holographic QCD models [95-98]. The scalar 
glueball ^ is associated with the local gauge-invariant QCD operator tr{G^yG^^'^) defined 
on the boundary spacetime, which has dimension Ac^ = 4. From the AdS/CFT dictionary, 
the scalar glueball has zero 5D mass, i.e. 5 = 0- 

We assume the glueball can be excited from the QCD vacuum described by the 
quenched dynamical holographic model in Section 2.1, and the 5D action for the scalar 
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n(U ) 


Soft-wall model 








HG = 430 


/iG = 600 


fiG = 1000 


1 


1216 


1697 


2828 


2 


1490 


2078 


3464 


3 


1720 


2400 


4000 


4 


1923 


2683 


4472 



Table 2. 0++ glueball in the original soft-wall model. The unit is in MeV. 



glueball ^{x,z) in the string frame takes the form as that in the original soft-wall model 
[95, 96] 

The only difference is that the metric structure in the original soft-wall model is AdSs, 
but in our dynamical soft-wall model the metric structure is selfconsistently solved from 
Section 2.1. 

The Equation of motion for has the form of 

- e-(3^-*)a,(e3^-*a,^„) = ml^^<^n- (2.24) 

After the transformation ^„ — )• e~'^^'^^^~^^'^m we get the schrodinger like equation of 
motion for the scalar glueball 

- < + V^'^n = (2.25) 
with the 5D effective schrodinger potential 

\ . (2.26) 

2.2.1 Glueball spectra in original soft-Yi^all model 

In the original soft-wall model for glueball [95, 96], the dilaton background takes the 
quadratic form <I> = i^qZ^ but the metric structure is still AdSs, one can easily derive 
the Regge spectra for scalar glueball: 

= 4/^g(^ + 1), n = 1, 2, • • • (2.27) 

which implying that the Regge slope for the scalar glueball is 4/2^., and the lightest glueball 
mass square is 8//^. In Table 2, we list some numerical results for the scalar glueball based 
on Eq. (2.27) with hg = 0.43,0.6, IGeV, respectively. 

From the lattice data for the scalar glueball as given in Table 1, one can read that 
the slope of the Regge spectra is around 4GeV^, which means /xg — IGeV. Prom Eq. 
(2.27), the lightest scalar glueball mass square in the soft-wall model should be around 
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n(0++) $ = fi^z^ 



2 




1 



3 



fiG = 900 HG = 1000 HG = 1100 
1434 1593 1752 
2356 2618 2880 
2980 3311 3642 
3489 3877 4264 



Table 3. 0++ glueball spect ra in the soft-wall model with positive quadratic dilaton background 
<f> = fi^z^ in unit of McV. 



gives //G — 0.5 GeV, then the slope for the Regge spectra will be around 1 GeV^, which is 
too small comparing with the lattice results 4 GeV^. 

In summary, by using the AdSs metric, the soft-wall model with the quadratic dilaton 
background field cannot accommodate both the lightest scalar glueball mass and the Regge 
slope. 

2.2.2 Glueball spectra in quenched dynamical hQCD model 

In the dynamical soft-wall (DSW) model, we firstly investigate the scalar glueball spectra 
with the positive quadratic dilaton background Eq.(2.17). 

Under the boundary condition ?^„(0) — )• and §f^(oo) — )• 0, we get the scalar glueball 
spectra as shown in Table 3. It is observed that with 0.9 GeV < hg < 1-1 GeV, the 
scalar glueball spectra in the dynamical soft-wall model with positive quadratic dilaton 
background can fit lattice results quite well. 

We would like to emphasize that the dynamical soft-wall model has the same parame- 
ters as the original soft-wall model, i.e. the AdSs radius L which is taken to be 1, and the 
quadratic coefficient of the dilaton background field /.ic. As we have shown in Sec. 2. 2.1, 
the original soft-wall model cannot accommodate both the ground state and the Regge 
slope. However, if one self-consistently solves the metric background under the dynamical 
dilaton field, it gives the correct ground state and at the same time gives the correct Regge 
slope. This is a surprise result! To explicitly see the difference, we show the scalar glueball 
spectra in the soft-wall model (blue dash-dotted line) and the dynamical soft-wall model 
(red solid line) in Fig. 1 for the case of fiQ = IGeV. 

It is observed from Fig. 1 that the glueball spectra in the dynamical soft-wall model 
is parallel to that in the soft-wall model, and the separation is about 5.8/i^. This indicates 
the ground sate of the scalar glueball has mass square around m|.^^^ = 2.5^^, and has 



mass around mcg^n=i = V^^^G■ From numerical results, we extract the Regge spectra in 

the dynamical soft-wall (DSW) model: 




^IL J. 

2 ~ 3GeV^ which 




2,DSW 



4^2^(„ _ 1) + 2.5^1, n = l,2,--- 



(2.28) 
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1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Figure 1. The 0++ glueball spectra for $ = /i^z^ with = IGeV in the soft-wall model (blue 
dash-dotted line) and the dynamical soft- wall model (red solid line) and compare with lattice data. 




Figure 2. The potentials and their difference. 



In order to understand the difference between the soft-wall model and the dynamical 
soft-wall model, we plot the effective schrodinger potentials Vc^ of the two models and their 
difference in Fig. 2. It is observed that the schrodinger potential (red solid line) in 
the dynamical soft-wall model has a lower minimum than that in the soft-wall model (blue 
dashed line), the difference is about 5.8//^, which is the same as the difference of the mass 
square in these two models, i.e. Vc^^dsw ~ V'^,sw = '^l' sw ~ ^% DSW ~ 5-8//q. 

If the dynamical soft-wall model takes the negative quadratic dilaton background $ = 
—^qz'^, the metric structure has the form of Eq.(2.18), and the scalar glueball spectra is 
shown in Fig. 3 with fic = 0.43,0.6, IGeV, respectively. 

It is observed that the negative quadratic dilaton background can also generate the 
Regge spectra. However, like the soft- wall model, the dynamical soft- wall model with 
negative dilaton cannot accommodate both the ground state and the Regge slope. 
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I: Latl 

t: Lal2 
:Liit3 

>: Lat4 
Lat5 



^^a = 1 GeV 



liG = 0.6aeV 



lAc, = 0.43 GeV 



1.0 1.5 



3.0 3.5 4.0 



Figure 3. 0++ glueball in tlie dynamical soft-wall with negative quadratic background $ 
with fiG = 0.43, 0.6, IGcV. The unit is in GeV and the dots are lattice data. 



2.3 Linear quark potential in quenched dynamical soft-wall model 

We follow the standard procedure [42, 99] to derive the static heavy quark potential VQQ^r) 
in the dynamical soft-M^all holographic model under the general metric background Eq. 
(2.4). In SU(N) gauge theory, the interaction potential for infinity massive heavy quark 
antiquark is calculated from the Wilson loop 



W[C] = ^TrPexp[z I Af.dx''] 



(2.29) 



where is the gauge field, the trace is over the fundamental representation, P stands 
for path ordering. C denotes a closed loop in space-time, which is a rectangle with one 
direction along the time direction of length T and the other space direction of length Rqq- 
The Wilson loop describes the creation of a QQ pair with distance R at time io = 
and the annihilation of this pair at time t = T. For T — t- oo, the expectation value of 



the Wilson loop behaves as {W{C)) oc e 



-TV, 



QQ. 



According to the holographic dictionary. 



the expectation value of the Wilson loop in four dimensions should be equal to the string 
partition function on the modified AdSs space, with the string world sheet ending on the 
contour C at the boundary of AdSs 



[W^'-iC]) 



^string i^] 



'Sng [C] 



where Sng is the classical world sheet Nambu-Goto action 

1 



Sng 



27ra„ 



(fr]^/Detx, 



abi 



(2.30) 



(2.31) 



with Up the 5D string tension which has dimension of GeV~^, and Xab is the induced 
worldsheet metric with o, b the two indices of the world sheet coordinates {rp , i]^). Without 
loss of generality, we can choose the if = t, r]^ = x, and the position of one quark is 



X 



and the other is x = — f^. Under the background (2.4), the Nambu-Goto 
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action Eq.(2.31) becomes 

TV 



Sng = ^ I dxe^^'Wl + z'^, (2.32) 



with the prime ' denotes the derivative with respective to x. 

Since there's no dependence on x, we can easily obtain the equation of motion: 



l\2 



Constant = e^-^^^^o), (2.33) 



for the minimum world-sheet surface configuration. 

Here the Rqq is dependent on zq which is the maximal value of z and z'{x = 0) = 0. 
For the configuration mentioned above and the given equation of motion, we impose the 
following boundary condtions z{x = 0) = zq,z{x = ±-f^) = 0. Following the standard 
procedure, one can derive the interquark distance Rqq as a function of zq 

R^oizo) = 2 r dz , ^ (2-34) 

V 

The heavy quark potential can be worked out from the Nambu-Goto string action: 

Vr^oizo) = ^ r dz^M^^, (2.35) 



with 9p = It is noticed that the integral in Eq.(2.35) in principle include a pole 
in the UV region (z — )• 0), which induces Vqq{z) — t- oo. The infinite energy should be 
extracted through certain regularization procedure. The divergence of Vqq{z) is related to 
the vacuum energy for two static quarks. Generally speaking, the vacuum energy of two 
static quarks will be different in various background. In our latter calculations, we will use 
the regularized ; which means the vacuum energy has been subtracted. A minimal 

subtracted result related to the background solution Eq.(2.15) is as following, 

Roo{zo) = 2zo / dv , ^ = ,y^°l - (2.37) 



The integrate kernel in Eq.(2.36) has a pole at = 1, and by expanding the integral 
kernel at u = 1 one has 

From Eqs. (2.36,2.37,2.38), we can find the necessary condition for the linear quark potential 
is that: There exists a point Zc, at which 

b'^{zc) 0, bs{zc) const, (2.39) 
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then the integral is dominated hy u = 1 region, one can obtain the string tension 




(2.40) 



Figure 4. The metric structure bs{z) = e"*=^^^ as functions of z corresponding to the solution 
Eq.(2.17) for $ = /i'qz'^ (red sohd line) and Eq.(2.18) for $ = —IIqz'^ (black dashed line), respec- 
tively. The blue dash-dotted line stands for the pure AdSs case, hg = IGeV has been taken for 
numerical calculation. 

Fig.4 shows the metric structure bs{z) as functions of z for the AdSs metric (blue dash- 
dotted line), and for the solutions Eq.(2.17) (red solid line) and Eq.(2.18) (black dashed line) 
of the quenched dynamical soft-wall model with positive and negative quadratic dilaton 
background fields $ = ±/x^z^ respectively. We can see that only for the case of positive 
dilaton background <I> = fJ-Qz'^, the metric solution Eq.(2.17) has a minimum point Zc- 
Therefore, the quark-antiquark potential should have a linear part for positive quadratic 
dilaton background $ = ij^'qZ^^ which can be seen explicitly from Fig. 5. While for the pure 
AdSs case as well as for the dynamical soft-wall model with negative dilaton background 
field $ = —fj.Qz'^, there doesn't exist a Zc where b'g{zc) — )• 0, and correspondingly the heavy 
quark potential does not show a linear behavior at large z. 

2.4 Short summary 

In this section, we have modeled the pure gluon system by using the quenched dynamical 
soft-wall model in the graviton-dilaton framework. The quadratic dilaton background field 
$ has dimension of 2, and is dual to the the dimension-4 gluon condensate. Comparing 
with the original soft-wall model with AdSs metric, here the metric background at IR is 
self-consistently deformed by the gluon condensate. 

It is found that the positive quadratic dilaton background can give the correct glueball 
spectra including the Regge slope and ground state, as well as the linear quark potential, 
and the negative quadratic dilaton background field can be safely excluded. In the following 
study, we will only focus on the case of positive quadratic dilaton background. 
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1.5 r 




Figure 5. The quenched quark potential result Vqq as functions of Rqq in the quenched dynamical 
soft-wall mode with positive $ = h'qZ^ (red solid line) and negative $ = —fi^z'^ (black dashed line) 
quadratic dilaton background fields, respectively. The blue dash-dotted line stands for the pure 
AdSs case, fic ~ IGeV and gp — 0.4 have been used for numerical calculations. 

3 Two flavor system: KKSS model and improved KKSS model 

We now turn to the the light flavor system with chiral symmetry SU{2)l x SU{2)pi. As 
we have mentioned in the Introduction, the current achievements of AdS/QCD models for 
hadron spectra are the hard- wall AdS/QCD model [21] and the soft-wall AdS/QCD or 
KKSS model [22] and its extended version [34-36]. In the hard- wall model [22], the chiral 
symmetry breaking can be realized by chiral condensation in the vacuum, however, the 
resulting mass spectra for the excited mesons behave as ~ n^, which is different from the 
linear Regge behavior ~ n. In order to generate the linear Regge behavior, the authors 
of Ref.[22] introduced a quadratic dilaton background, one can obtain a desired mass 
spectra for the excited vector mesons, while the chiral symmetry breaking phenomenon 
cannot consistently be realized. In the following, we firstly give a brief introduction on the 
KKSS model and review the meson spectra in this model. 

3.1 The KKSS model 

The KKSS model [22] have two background fields: the positive quadratic dilaton back- 
ground <I> = /i^z^ and the metric background qmn- Note, in the following, we will use /U 
instead of to distinguish from the pure gluon system. The background geometry is not 
dynamically generated but assumed to be AdSs space with the metric structure 

ds'^ = QMNdx^^dx^ = -2 {'nf,^dx''dx'' + dz^) , (3.1) 

which gives As{z) = — log(z/L). 

The mesons are described by 5D fields propagating on the background with the action 
given by 

Skkss = - j d^xe-''^^^^sTr{\DX\'' + mlX^ + ^{Fl + Fl)), (3.2) 



-14- 



with 55 = VItP' jNc- The scalar field X is dual to the dimension-3 qq operator, and rax 
is the 5D scalar mass. According to AdS/CFT dictionary, the dimension-3 scalar has 5D 
mass = —3. The field X[z) is actually a complex field to incorporate the scalar S and 
the pseudoscalar P fields, 

X"/'(^)=/'^ + 5^r/^e^i^"*", (3.3) 



where a, /3 are in the isospin space, a = 1,2,3 are the SU{2) generator index. The scalar 
field takes a nonzero vacuum expectation value (VEV) xi^)i which is expected to realize 
the chiral symmetry breaking. 

The Gauge fields Lm and Rj^j model the SU(2)lx SU(2)^ global chiral symmetry of 
QCD for two flavors of quarks, which are defined as 

where L^'^ = L^^°-t°- and Tr[t°'t^] = (5"^/2. The covariant derivative becomes 

jjMj^ = d^^X - iL^X + iXR^^. (3.5) 

To describe the vector and axial-vector fields, we simply transform the L and R gauge 
fields into the vector (V) and axial-vector (A) fields with L*^ = + and i?*^ = 
-A^, one can have F ^ + F I = 2 (F^ + Fj) , with 



j,M IV ^ QM ylV _ QIV yM _ _ ^yM ^ylV^^ (g^g) 

v2 

pMN ^ qMj^N _ qN^M _ J^[^A/^ A^]. (3.7) 

v2 

In terms of the vector V and axial- vector A fields, the KKSS action Eq.(3.2) can be 
rewritten as 



Skkss = - j d^x^e-*(^)Tr 



|Z)X|2+m^|Xp + -L(F2+Fl) 



(3.8) 



where the covariant derivative now becomes 

D^^X = d^^X - i[V^'^,X] - i{A^',X}. (3.9) 
3.2 Degeneration of chiral partners in KKSS model 

The scalar field takes a nonzero vacuum expectation value (VEV) x(z), which is expected 
to realize the chiral symmetry breaking as in the hard wall model. We will show in the 
following that the chiral symmetry breaking is not realized in the soft-wall model or KKSS 
model, and we will analyze the reason. 

Scalar vacuum expectation value 
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The equation of motion for the scalar vacuum expectation value (VEV) xi^) defined 
in Eq.(3.3) can be deduced and takes the following form, 

x" + (3^1 - $')x' - mj^e^^^x = 0. (3.10) 

In the hard wall model, <&' = 0, the scalar VEV has the exact solution 

X{z) = ciz + C2Z^ = irigZ + az^, (3.11) 

where we have identified niq = c\ and a =< qq >= C2. In the softwall model, ^{z) = fi'^z'^ 
and = 2//^z, and the general solution of Eq.(3.10) has the form of 



X{z) = c^Gl'' { -z'\ ,\ ] +c,e^z^ f/o ('^\ +h ('4 ] ) • (3-12) 



2' 2 



with In{z) the modified Bessel function of the first kind,G^" 





ai, . . 


.,ap\ 






bi,. 





the MeijerG 



function(For details, please refer to Mathematica 8.0). 

To be more instructive, we can extract the large z behavior of x from the equation of 
motion : assuming x" << Xi^hen z >> 1, we have —2iJ?zxi + 3x/2^ = 0, and assuming 
x' >> Xi^hsn z » 1, we have x " ~ 2^^zx' = 0. And then we could get the IR behaviors 
of the two independent solution: xi ~^ e^^^"^ /{pLz) and X2 — ^ e^-^/i'^t^^^'^) — ). i. The first one 
would make the spectra of oi nonlinear, so in order to produce linear ai spectra, ci, C2 in 
Eq.(3.10) are not independent. Requiring x oc X2(^), when z » 1, the small z expansion 
of X would be: 

X(^) = c(^z + ^3^3 1^-1 + 7s + ttH + iog(^^)^ ) + 0{z') (3.13) 

= c{fiz + /i^ 2:^(0.095 + log(/xz))) (3.14) 

with = 0.577 the Euler's constant and iIj{z) the digmma function with ^(— ^) = 0.036. 
We would take c = mq/fj,, so a can be read as o" = 0.095mq/i^. The experimental data for 
vector, axialvector, scalar and pseudoscalar are shown in Table 4. To fit the Regge slope 
of vector meson p, we have to choose p = 0.43. Then even we take niq = 9MeV, a is only 
(54 MeV)^, which is too small comparing with the experienced value (250 MeV)^. This 
problem was pointed in the original paper and the authors also mentioned to add quartic 
terms \X\^ to cure it. 

Meson spectra 

In the following, we show the meson spectra in the KKSS model. The equations of 
motion of the vector, axial- vector, scalar and pseudo-scalar mesons take the form of: 

- P'n + VpPn = mlpn, (3.15) 

-a„ + Vaan = mlan, (3.16) 
-s'^ + VsSn = mlsn, (3.17) 

+ V^ifn = gle^'xi^n - e^'x^n)- (3.18) 
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n 


p (MeV) 


ai (MeV) 


/o (MeV) 


vr (MeV) 


1 


775 ± 1 


1230 ± 40 




140 


2 


1282 ± 37 


1647 ± 22 


980 ± 10 


1300 ± 100 


3 


1465 ± 25 


1930^^5] 


1350 ± 150 


1816 ± 14 


4 


1720 ± 20 


2096 ± 122 


1505 ± 6 


2070 


5 


1909 ± 30 


2270l|^ 


1724 ± 7 


2360 


6 


2149 ± 17 




1992 ± 16 




7 


2265 ± 40 




2103 ±8 




8 






2314 ± 25 





Table 4. The experimental data for meson mass from PDG [31]. The data selection scenario used 
here is the same as in Ref.[35], which shows the chiral symmetry breaking maintains in the highly 
excited states of chiral partners. 



n 


p (MeV) 


ai (MeV) 


/o (MeV) 


TT (MeV) 


1 


860 


860 


1053 


1054 


2 


1216 


1217 


1360 


1360 


3 


1489 


1490 


1609 


1609 


4 


1720 


1720 


1824 


1824 


5 


1923 


1923 


2017 


2017 


6 


2107 


2107 


2193 


2192 


7 


2275 


2275 


2355 


2355 



Table 5. The mass spectra for vector mesons p, axial vector mesons ai, scalar mesons /o and 
pseudoscalar mesons tt in the KKSS model with = 9 MeV,/U = 430 MeV, which gives a = 
(54 MeV)3). 



with schrodinger like potentials 



^P = ^^V^+ ^ ^ , (3.19) 

Va = + ^^^-f^" + gle'^^x\ (3.20) 

y. = M^+(M.^_^2^e-^ (3.21) 



^ _ 3^: - ^" + 2x"/x - 2x'Vx' , (3^1 -<^' + 2x/x? ^ (3^32) 



= ^^^^ + (3.23) 

Since g^e^^'x'^ ~^ when z — t- 00, we can expect the spectra of the chiral partners, i.e. the 
vector and axial vector as well as the scalar and pseudoscalar mesons would be degenerate 



-17- 




4 



-2 



Figure 6. Meson spectra in the KKSS model with vriq = 9 McV, \Ji = 430 MeV comparing with 
experimental data in Table 4. 



in the large n region. 

The meson spectra in the KKSS model is shown in Table 5 and Fig. 6. In order to realize 
the linear Regge behavior, we have used parameters ruq = 9 MeV, ^ = 430 MeV. However, 
this gives a small chiral condensate a = (54MeV)^, which leads to the degeneration of 
chiral partners, i.e. the scalar spectra overlaps with the pseudoscalar spectra, and the 
vector spectra overlaps with the axial-vector spectra. On the other hand, in order to 
realize the chiral symmetry breaking in the KKSS model, i.e. the separation of the spectra 
of the chiral partners, as shown in [33], one cannot get the linear Regge behavior for the 
axial- vector meson. 

3.3 Improved KKSS model with quartic interaction term 

As we have shown above that the KKSS model cannot accommodate chiral symmetry 
breaking and linear confinement. Refs. [34, 35] introduced a quartic interaction term kX'^ 
in the bulk scalar potential to improve the situation. Nevertheless, such a term was shown 
in Ref. [35] by Gherghetta-Kapusta-Kelley to result in a negative mass for the lowest lying 
scalar meson state. 

The meson spectra in the Gherghetta-Kapusta-Kelley (GKK) model is shown in Fig. 7, 
where the parameters are chosen: rUq = 9.75MeV,o- = (204.5 MeV)^,/^ = 430 MeV (equiv- 
alent to A = 0.183Gey2 in their notation). The lowest scalar meson has mass square 

n=i ~ — 0.559GeV^, which shows the instability in the scalar sector. 

3.4 Improved KKSS model with deformed warp factor 

In Ref. [36], Sui-Wu-Xie-Yang (SWXY) introduced a deformed warp factor in the KKSS 
and GKK model and the metric structure takes the form of 



which can cure the instability of the scalar potential and produce meson spectra in good 
agreement with experimental data. Even though the authors in Ref. [36] didn't calculate 



bs{z) 




(3.24) 
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Figure 7. Meson spectra in tlie GKK model witli niq = 9.75 MeV, a = (204.5 MeVf,fi = 430 MeV. 



the heavy quark potential, but from our criteria for the linear quark potential Eq. (2.39), 
the geometric factor Eq.(3.24) in the SWXY model can produce a linear potential. 
In this model, the authors grouped their settings into "case-a" and "case-b". 

1) In "case-a", the large z behavior of the scalar (see Eq.(9) in their paper) is xi^ ~^ 
oo) = j{fiz), since their metric warp factor is like As{z — > oo) = constant, the difference 
between the effective potential in vector sector and axial-vector sector takes the limit of 
d'ix^^'^^' oc z^, we can see the Regge slope for the vector spectra is different from that for 
the axial- vector mesons. 

2) In "case-b", x{z — )■ oo) = the difference between the effective potential in 
vector sector and axial-vector sector takes the limit of glx^e^^" oc z, and the vector and 
axial-vector spectra would approach each other at high excitations. 
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( case-a ) 



(case-b ) 



Figure 8. Meson spectra in the SWXY model comparing with experimental data in Table 4. The 
parameters ^ = 445MeV (their fid), and niq = 4.98MeV, a = (255MeV)^ are used for case-a, and 
rriq = 4.25MeV, a = (268MeV)3 are used for casc-b. 

The meson spectra for the SWXY model are shown in Fig.8(we only take the Model 
III in their original paper as an example), where for case-a(their Ilia) they have used 
parameters as /x = 445MeV,??T,q = 4.98MeV,iT = (255MeV)^, and for case-b(their Illb), 
they have use parameters as /u = 445MeV,mg = 4.25MeV,(T = (268MeV)^. It is found 
that for case-a, the Regge slopes for the scalar and vector meson spectra are the same, and 
the Regge slopes for the pseudo-scalar and axial-vector meson spectra are the same, while 
the slopes of the two groups are different. For case-b, all spectra are degenerate. 

It is worthy of mentioning that the meson spectra are compared with experimental data 
taken in Table 4, which are the same as in Ref. [35], and different from the experimental 
data taken in Ref. [36]. As for which data should be taken properly, and whether there 
should be chiral symmetry restoration at high excitation states [33, 100, 101], we leave 
them as open questions. 

4 Two flavor system: the dynamical soft-wall model 

A successful holographic QCD model should describe chiral symmetry breaking, and at the 
same time should describe both the Regge trajectories of hadron spectra and linear quark 
potential, two aspects in the manifestation of color confinement. Thus how to naturally 
incorporate all these important features into a single AdS/QCD model and obtain the 
consistent mass spectra remains a challenging and interesting task. In this section, we 
provide a fully dynamical soft-wall holographic QCD model formulated in the graviton- 
dilaton-scalar system, which can incorporate chiral symmetry breaking, Regge spectra as 
well as linear quark potential. 



-20- 



4.1 Dynamical soft- wall model: the graviton-dilaton-scalar system 

As we have shown in Sec. II that the pure gluodynamics can be described very well by 
the quenched dynamical soft-wall model formulated in the graviton-dilaton system. The 
quadratic dilaton background related to the gluon condensate in the vacuum can produce 
the linear confinement, including linear Regge spectra and the linear heavy quark potential. 
We now add light flavors in terms of meson flelds on the gluodynamical background. The 
total 5D action for the graviton-dilaton-scalar system takes the following form: 

S = Sg + ^^Skkss, (4.1) 

with 

Sg = / d'x^e-^''{R + ^dM^d^^ - Vci^)), (4.2) 

Skkss = - j d'x^se~''Tri\DX\^ + VxiX+X, $) + ^(^1 + ^I))- (4-3) 

It is noticed that Sg is the 5D action for gluons in terms of dilaton field ^> and takes the 
same form as Eq.(2.3), here we have assumed the action is in the string frame. Skkss is 
the 5D action for mesons propagating on the dilaton background and takes the same form 
as the general KKSS action Eq.(3.2). Vg'(^) and VxiX~^ X,^) are potentials for dilaton 
field and scalar field, respectively. It is noticed that the scalar field might mix with the 
gluon fields, therefore we have chosen a general form for the scalar potential Vx{X^ X, $). 

4.2 Background flelds in the vacuum with chiral and gluon condensate 

In the vacuum, it is assumed that there are both gluon condensate and chiral condensate. 
The square of the dilaton background field is supposed to be dual to the gauge invariant 
dimension-4 gluon condensate < >. In order to distinguish the dilaton background from 
the pure gluon system, we take the positive quadratic dilaton background with a different 
coefficient 

$ = n^z^. (4.4) 

The scalar field X{z) is a complex field as shown in Eq.(3.3) and it is expected that the 
scalar field takes a nonzero vacuum expectation value (VEV) x{^)- 
It's easy to get the 5D action for the vacuum background: 

Svac = SG,vac + -irr S K K S S ,vac, (4.5) 



with 



SG,vac = / d'x^se-^'' {R + AdM^d^^ - Vg{^)) (4.6) 

SKKSS,vac = - j <fx^se-'^{\dMXd^'x + ^c(x, ^)) (4-7) 
where we have defined Vc = TriVx) For further convenience, we define 

^<^.x-^> ^c,xx--^- (4-8) 



- 21 - 



By Redefinition: La ;y — )■ L^Vc — )• Vc, ^^^3^^ ^ — )■ A we can set all the fields and constants 
to be dimensionless, and the vacuum action takes the form of 

Syac = J d'x ^sie-^^lR + ^dM^d^'^ - Vein 

-Xe"''{^dMXd''x + Vb(x, <&))} (4.9) 

After the frame transformation ^f^^ = gfjj^fe^^ , the action Syac in the Einstein frame 
takes the following form 

Svac = / d^'x ^{[Re - pM^d''^ - vim 

-Xe^i^dMXd^'x + el'^Vcix, ^))}. (4.10) 
The Einstein equation and field equations in the Einstein frame have the expression of 

Emn + \giiN [^mO''^ + Vi{<^) + A(ie*aa9'x + el^Vcix, ^))) (4.11) 

-pM^dN^ - ^e^'dMxdNX = 0, (4.12) 
^ dMiVa^d''^) - xle''dMXd''x - (Vg{^) + Aei*Vb(x, f )) = 0, (4.13) 



3^5^ ^"-^ ' 2 

X^dMiVg^e^d^x) - 5x (^g(^) + Xel^Vcix, ^)) = 0. (4.14) 



/9E 

We can derive the three coupled field equations in the string frame as 



< + ^? + V-lA>'-^eV = 0, (4.15) 



$ ■ + (3^, - 2$ )<!>■ - ^e^x' - ^e2^-t*a$ (Vg{'^>) + Aei*Vfc(x, f ) j = 0, (4.16) 

X + (3^1 - ^')x' - e^^'Vc-yix, ^) = 0. (4.17) 

If we know the form of the dilaton field $ and the scalar field Xj then the metric Ag, the 
dilaton potential Vg{^) and the scalar potential Vc{Xi ^) should be self-consistently solved 
from the above three coupled equations. 

4.3 Chiral symmetry breaking and linear confinement 

In the case of positive quadratic dilaton background ^{z) = /i^z^, we will constrain the 
form of scalar VEV from the linear confinement. 

The UV asymptotic form of xi^) 

As proposed in the KKSS model, at the ultraviolet (UV) region, the scalar field takes 
the following asymptotic form, 

xiz)'-^m,Cz + ^z^ (4.18) 
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where niq is the current quark mass, and a is the quark antiquark condensate, and C is a 
normahzation constant and is fixed as C = ^^a^^. with Nc = 3, Nf = 2. 

The IR asymptotic form of x{^) constrained from linear quark potential 

The hnear behavior of quark-antiquark static potential in the heavy quark mass hmit 
mq — >• oo can describe the permanent confinement property of QCD. In Sec. 2.3, we have 
derived the heavy quark potential under the general metric background As, and we have 
observed that for the metric structure hs = e^" , if there exists a point Zc where hg{zc) — )• 0, 
then one can extract the string tension as of the linear potential as 

a. = |iM^^-f6^(..). (4.19) 

Therefore, the necessary condition for the linear part of the Q — Q potential is that there 
exists one point Zc or one region, where bs{z) — t- 0, z — )• while bs{z) keeps finite. For 
simplicity, we can take the following constraint on the metric structure at IR: 

A'siz) 0,As{z) Const. (4.20) 

Under the condition of Eq.(4.20), the equation of motion Eq.(4.15) in the IR takes the 
following simple form: 

- = 0, (4.21) 

6 b 

which provides a relation between the chiral condensate and low energy gluodynamics at 
IR. The asymptotic form of xi^) &t IR can be solved as: 



Xiz) ysTA/ze-*/^. (4.22) 

The full form of xi^) 

To match the asymptotic forms both at UV and IR in Eqs.(4.18) and (4.22), x can be 
parameterized as 

xiz) = 787A/ie^*/2(l + cie"* + cae-^*), (4.23) 

with ci = -2 + '-^^ + and C2 = 1 - '-^^ - Solving Eq.(4.23), we can 

obtain the full expression for the scalar VEV, which takes the following form: 

X{z) = ^5^\/^(5^Erf (-8^2(^3 + + b^^Xfi' 

+3(\/5Erf ^yi^i^^ (4\/2C/i^ - 6\^a - ^C^VXfi^mq 
+20^/2C^3Erfj^^^)), (4.24) 

where gl = Air^^ and = 
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With the input of ^{z) = ^^z^ and x{z) given in Eq.(4.24), one can solve the metric 
As or hs from the equation of motion Eq. (4.15). By taking G5 = ^^,mq = 5 MeV,cJ = 
(240 MeV)^, we show the numerical results for the scalar VEV x{z) ^'^^ the solved metric 
structure hs{z) in Fig.9. 




( X ) {bs) 



Figure 9. Scalar VEV x(z) and solved metric structure hs as functions of z with G5 = niq = 
5 MeV,cr = (218 MeV)3. 

The heavy quark potential under the solved metric structure is also shown in Fig. 10 
by the solid line and comparing with the Cornell potential yComeii^j^^ = _|_ (j^ji ^ 
with K 0.48, as ~ 0.183GeV^ and Vq = — 0.25GeV. It is observed that the heavy quark 
potential produced in our model including the back-reaction from light flavor dynamics 
agree well with the Cornell potential. 




Figure 10. Heavy quark potential Vqq as a function of Rqq from our model (solid line) with 
gp = 2.3 and /i = 430MeV, a = (218MeV)^ compared with the Cornell potential (dot-dashed line). 



4.4 Meson spectra in the graviton-dilaton-scalar system 

With the dilaton background field ^{z) = fi'^z'^ and the scalar background field xi^) given 
in Eq.(4.24), we have solved the metric As or bg from the equation of motion Eq. (4.15). 
Now we are ready to derive the meson spectra in the dynamical soft-wall model. 

4.4.1 Scalar spectra 

The action of the scalar perturbation is 

Ss = J d'xe-^'^sid.Sd'S + d^Sd'^S + Fc,xx(x, ^)S^), (4-25) 

and the equation of motion for the scalar perturbation after doing the KK modes expansion 
is 

- e-(3^-*)a,(e3^-*9,5„) + e2^^yc,xx(x, ^)sn = mls^. (4.26) 
By doing the transformation Sn — )• SnC"^^"^^^*^^^, one can get the schrodinger like equation 

- Sn^Vs[,z)sn = m\sn (4.27) 

with the schrodinger potential 

3<-ci." , (3a;-c^')2 



ys{z) 



+ 



+ e'^°^c,xx(x,$). 



2 4 
Assuming the scalar potential can be separated into 

Vb(x,$) = e^(*V,(x), 



(4.28) 



(4.29) 
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rUg (MeV) 


a^/^ (MeV) 


Mod I 


0.75 


5 


218 


Mod II 


0.75 


5 


240 



Table 6. Two sets of parameters. 



we have 



X 



(4.30) 



By using Eq.(4.17), we can have 



^ ^c,xx ~ ^ I 

X 

III tl 

X X 
-(2^1 + /,$$')(3^1 -<!>'). (4.31) 

In our solution, x[A oc e~2*j so the leading term of e^^'^Cxx ^ ^("^ ~ 

/,*) + /,$*^'^ = (f + ^4^)*^'^- The Regge behavior of the spectral is determined by the 
leading IR behavior, the coefficient before is proportional to the Regge slope. In order 
to be consistent with the experimental data, we need | + — )• in the IR region, so 
the leading term of = —\ and /($) — )■ — Then we examine the IR behavior of the 
EOM of X, 

- ye-t - $'e-t oc e^(*Ve,^(x) (4.32) 

Note that when z — t- oo, we have $' oc 2; oc and x ~^ const, Vc^y.{x) — ?• const, then we 
can know that when $ — )• 00, /(<!*) = + ^log<I>. If we hope in the small <I> region, 
g/(*) = 1^ a simplest choice is 

/(*) = -f + i5£^^. (4.33) 

which leads to the coupling between dilaton background field and the scalar field at leading 
order taking the form of 

^c(x,<J>)~x'«''- (4.34) 

For our numerical calculations, we take two sets of parameters in Table 6. The param- 
eters in Mod I has a smaller chiral condensate, which gives a smaller pion decay constant. 
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n /o Exp (MeV) Mod I (MeV) Mod II (MeV) 



1 


cc;n+250 


460 


231 


2 


980 ± 10 


1069 


1064 


3 


1350 ± 150 


1381 


1395 


4 


1505 ± 6 


1629 


1632 


5 


1724 ± 7 


1845 


1846 


6 


1992 ± 16 


2038 


2039 


7 


2103 ± 8 


221 


2215 


8 


2314 ± 25 


2376 


2376 



Table 7. The experimental and predicted mass spectra for scalar mesons /q. 




( Mod I ) 



( Mod II ) 



Figure 11. Scalar meson spectra nij^ „ as functions of n for Mod I and II defined in Table 6. 

this will be discussed in Sec. 5. With two sets of parameters, the scalar spectra are shown 
in Table 7 and Fig. 11, and the corresponding wave-functions are shown in Fig. 12. 

It is observed that in our graviton-dilaton-scalar system, the lowest scalar state has 
a positive mass, and the higher excitations behave a Regge line which agrees well with 
experimental data. 



4.4.2 Pesudo-Scalar Sector 

The terms of quadratic order in vr and ip {A^jj_ = d^^p) is 



C(2) 



N 



f 



at 



(4.35) 
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Figure 12. The scalar wave function ipfg^niz) as function of z witli Gs = ^^^niq = 5MeV, ct = 
(218MeV)3. 



The equations of motion for the pesudoscalar vr coupled with (f take the form of 

- e-(3^=-*)9,(e3^^-*x'5,)^ + qWtt - ^) = 0, (4.36) 
-e-(^-*)a,(e^-*a,)(^-55Ve'^^(vr-<^)= 0, (4.37) 

which can be written in the following form 

- ip'^ + V^^n = gh'^'xi'^n-e'^^x^n), (4.38) 
with the effective schrodinger potentials 

, . _ 3A':-<^" + 2x"/x-2xVx' ^ jSA',-^' + 2x'/x? 

= ^^^^ + ^""^ ■ (4.39) 

With two sets of parameters in Table 6, the pseudoscalar spectra vr are shown in 
Table 8 and Fig. 13, and the corresponding wave-functions are shown in Fig. 14. It is 
observed that in our graviton-dilaton-scalar system, the lowest pseudoscalar state has a 
mass around 140MeV, which can be regarded as the Nambu-Goldstone bosons due to the 
chiral symmetry breaking. The higher excitations behave a Regge line which agrees well 
with experimental data. 
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n vr Exp (MeV) Mod I (MeV) Mod II (MeV) 



1 


140 


137 


139.4 


2 


1300 ± 100 


1498 


1600 


3 


1816 ± 14 


1865 


1897 


4 


2070 


2099 


2116 


5 


2360 


2287 


2299 



Table 8. The experimental and predicted mass spectra for pseudoscalar mesons tt. 



ffl-,(GeV-) »r.(GeV-) 




< i n < i 

2 4 6 8 2 4 6 8 



-2^ -2^ 

( Mod I ) ( Mod II ) 



Figure 13. Pseudo-scalar spectra „ as function of n for Mod I and II defined in Table 6. 



4.4.3 Vector sector 

In the vector sector, the terms of quadratic order in V-^ are 

4'^ = J d'xe-Hlid^V^^d^V^'^ + d.V^^d^V^''), (4.40) 

The equations of motion of the vector mesons take the form of 

- P'n + ^vPn = mlpn: (4.41) 

with the schrodinger like potential 

A' - (A' - 
K = + ^ ' ^ ' . (4.42) 

With two sets of parameters in Table 6, vector spectra are shown in Table 9 and Fig. 
15, and the corresponding wave-functions are shown in Fig. 16. It is observed that in our 
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Figure 14. The pseudoscalar wave function '4'T!,n{z) as function of z witli = ^^,mq = 
5MeV,cr = (218MeV)3. 



n p exp. (MeV) Mod I. (MeV) Mod II. (MeV) 



1 


775.5 ± 1 


750 


771 


2 


1282 ± 37 


1139 


1143 


3 


1465 ± 25 


1428 


1670 


4 


1720 ± 20 


1667 


1878 


5 


1909 ± 30 


1876 


2065 


6 


2149 ± 17 


2063 


2237 


7 


2265 ± 40 


2235 


2396 



Table 9. Tlic experimental and predicted mass spectra for vector mesons p. 



graviton-dilaton-scalar system, the lowest vector state has a mass around 770MeV, and 
the higher excitations behave a Regge line which agrees well with experimental data. 

4.4.4 Axial Vector Sector 

The terms of quadratic order in A-^ is 

= -^^2^ / d'xe-Hlid^A^d^A^f^ + d.A^d^A^^ + ^A^A^^^). (4.43) 
The equations of motion of the axial-vector mesons take the form of: 

- a'^ + Vattn = mlan, (4.44) 
with the schrodinger potential for the axial vector as 

Va = + + gy^^x'. (4.45) 
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As we can see that the difference between the schrodinger potentials for the axial vector 
Eq.(4.45) and vector Eq.(4.42) is only the extra term gle^^'x^ in Eq.(4.45). In the KKSS 
model, gle^^'x^ — ?• when z — )• cxd, therefore there is no splitting between vector and axial 
vector. Here in the graviton-dilaton-scalar system, g^e^^^y^ — )• constant when z — )• oo. 
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n ai Exp (MeV) Mod I (MeV) Mod II (MeV) 



1 


1230 ± 40 


1219 


1316 


2 


1647 ± 22 


1691 


1735 


3 


1930lf[3 


1948 


1969 


4 


2096 ± 122 


2149 


2163 


5 


22701^^ 


2326 


2336 



Table 10. The experimental and predicted mass spectra for axial vector mesons ai. 



iii-„,(GeV-) »r.,(GeV-) 




— ^ — , — , — \ — , — , — , — I — , — , — , — I — , — , — , — \ — , — ,„ ()_^ — , — , — \ — , — , — , — I — , — , — , — I — , — , — , — \ — , — ,j 

2 4 6 8 2 4 6 8 

-2L -2 . 

( Mod I ) ( Mod II ) 

Figure 17. m^^ „ as functions of n for Mod I and II defined in Table 6. 

which naturally induces the separation of the chiral partners. 

With two sets of parameters in Table 6, vector spectra are shown in Table 10 and 
Fig. 17, and the corresponding wave-functions are shown in Fig. 18. It is observed that 
in our graviton-dilaton-scalar system, the lowest axial vector state has a mass around the 
experimental value 1230Mey, and the higher excitations behave a Regge line which agrees 
well with experimental data. 

4.5 Short summary 

We summarize the meson spectra for scalar, pseudoscalar, vector and axial-vector in Fig. 
19 for two sets of parameters. It is found that for both cases, the results are in very well 
agreement with experimental data. The ground state has the order of m,r < "^/o < '^p < 
, and the Regge slopes for scalar, pseudoscalar, vector and axial- vector meson at high 
excitations take the same value of 4//^. 
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_2l -2l 

( Mod I ) ( Mod II ) 

Figure 19. All meson spectra in the dynamical soft-wall model with two sets of parameters in 
Table 6 comparing with experimental data. 

5 Decay constants, pion form factor and vector couplings 

So far, within the above holography model we have studied one aspect of the static hadronic 
properties, the resonance masses. To further confirm this model, we have to check that 
whether it can reproduce reasonable behavior of other static and dynamic properties of 
hadronic physics, such as decay constants, vector couplings and form factors etc. 
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As in Ref.[21] (for more details see Ref. [102]), by studying the current-current corre- 
lation function and rewriting it as a summation over the normalizable wave functions, we 
can extract the decay constants fn, Fp^^, Fa^^n as following, 

/- =-^e^=-*9,A(0,z)U^o, (5.1) 

Pi = ^(e^-*a.K(z)U^o)^ (5.2) 



Fl,n = ^(e^-*a,A„(z)|,^o)'. (5.3) 



Where ^4(0, z), Vn{z),An{z) is the solution of equations 

(_e-(^^-'f)5,(e^^-*a,) + gb'e^^^x^)A{id, z) = 0, (5.4) 

(_e-(^-*)5,(e^-*a,) - ml^)Vn{z) = 0, (5.5) 

(_e-(^-*)a,(e^-*9,) + gy^^x' - ml^^Mz) = 0, (5.6) 

with the boundary conditionA(0, 0) = 1,9^^(0,00) = 0, K(0) = 0,(9^K(oo) = 0,A„(0) = 
0,dzAn{oo) = and normalized as / dze^''~'^VmVn = J dze^'~^ AmAn = ^mn- 

We can also extract the pion form factor from the three point correlator as [103-106] 

PMQ") = 4^ f dze^^-''V{q',z){id.^f+glx'e'^^{7r-^)'}, (5.7) 

whereQ^ = — g^, and V{q'^ , z),TT{z),ip{z) is the solution of 

(_e-(^-*)5,(e^^-*a,) + q')V{q^ z) = 0, (5.8) 

_e-(3^-'^)5,(e3^-V5.)7r - „x'(vr - V') = 0, (5.9) 

_e-{A.-</.)5^(gA.-0^^)^ - <7ix'e2^^(7r - (^) = 0, (5.10) 

with the boundary condition V{q^,Q) = I, dzV{q^ ,oo) = 0, 7r(0) = Q,dz'n:{oo) = 0, (/^(O) = 
0, 99(00) = and normalized as 



N 



[ dze^^-mdz^f + glx'e'^'{^-^f} = l. (5.11) 

To make sure that -F7r(0) = 1 or equivalently we can write 

Jdze^^-^V{q^z){idz^r+glx'e'^^{^-V?} 
^ fdze^^-^{{dz^)^+glx^e^^<^-^?} ' 

Decomposing F^^ as in [103-106], we reach 

FAQ') = E^f^^ (5.13) 
Q'^ + mi, 

with 

J dze^^-'^Vn{{dzV? + glx^e'^^jT^ - ^f] 

9n.. 55 J-rf^eA.-$|(a^^)2+^2^2e2A.(^_^)2} ' ^5.14] 
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exp. (MeV) Mod I. (MeV) Mod II. (MeV) 

92.4 ±0.35 75 83^6 
F^^ 346.2 ±1.4 277 282 
FaC^ 433 ± 13 425 452 
gp^^ 6.03 ± 0.07 3^59 3.14 

Table 11. Decay constant in dynamical soft- wall model with two sets of parameters in Table 6. 





exp. (MeV) 


Mod II. 


KKSS 


hard-wall 


mod-soft 


SWXY 


GKK 


u 


92.4 ±0.35 


83.6 


87.0 


92.1 


88.0 


92.4 


92.4 




346.2 ± 1.4 


282 


261 


329 


325 








433 ± 13 


452 


558 


463 


474 






9p7V7V 


6.03 ± 0.07 


3.14 


3.33 


4.48 


4.63 


3.51 


2.89 



Table 12. Table for decay constants and couplings from other models, the results in hard- 
wall, KKSS model stands, mod-soft are taken from [105, 106]. 



We would denote gp^TT = go-Kn, i-e. putting the p meson ground state wave function Vp = Vq 
in the above equation. 

The numerical results for the decay constant are shown in Table 11 and are compared 
with other models in Table 12. 

We can see that both Mod I and Mod II give a smaller value of /tt than the experimental 
data, however we find that m^, <t, 771,^, /tt in Mod I, II obeys the Gell-Mann-Oakes-Renner 

2mqCJ = m'^f^ with — 1 of order 0.1% within numerical errors and it checks the 

1/2 1/2 

consistency of our model. The oi decay constants Fp ,Fal prediction is closer to the 

experimental data than other models, and both the prediction of the pion form factor are 

1/2 

better than the original soft- wall model in Ref.[105, 106], while Fp ,gpTT-K are a little too 
small. 

The pion form factor is shown in Fig. 20, it is found that for parameters used in Mod 
I with a smaller chiral condensate, which corresponds to a smaller pion decay constant, the 
pion form factor matches the experimental better. 

6 Discussion and summary 

In this work, we construct a quenched dynamical holographic QCD (hQCD) model in the 
graviton-dilaton framework for the pure gluon system, and develop a dynamical hQCD 
model for the two flavor system in the graviton-dilaton-scalar framework by adding light 
flavors on the gluodynamical background. The dimension-2 dilaton background field takes 
the same form as that in the soft-wall model, i.e. $ = jx^z^, and its square <I>^ can 
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_j QiiQ^yi 0.0 1 e'/GeV- 

5 1 2 3 4 5 



( Mod I ) ( Mod II ) 



Figure 20. Ft^{Q'^) as function of for Mod I and II defined in Table 6 and compared with 
experimental data. The blue dashed lines are the prediction in our model, and the green dotted line 
is the original soft-wall model results taken from Ref.[105, 106]. 

be identified as the dual field to the gauge invariant operator of the dimension-4 gluon 
condensate < TrG^ >. The gluon condensate in the QCD vacuum breaks conformality 
and naturally induces a deformed warp factor, which can be self-consistently solved in the 
graviton-dilaton and graviton-dilaton-scalar framework. 

In the quenched dynamical model, without introducing extra parameters but just self- 
consistently solving the deformed metric induced by the dilaton background field, we find 
that the scalar glueball spectra is in very well agreement with lattice data, while the soft- 
wall model with AdSs metric cannot accommodate both the ground state and the Regge 
slope for the scalar glueball spectra. We also give a necessary condition for the existence of 
linear quark potential from the metric structure, and we show that in the graviton-dilaton 
framework, a negative quadratic dilaton background field cannot produce the linear quark 
potential. 

For two flavor system in the graviton-dilaton-scalar framework, the deformed metric 
is self-consistently solved by considering both the chiral condensate and gluon condensate 
in the vacuum, which are responsible for the chiral symmetry breaking and linear confine- 
ment, respectively. It is found that the mixing between the chiral condensate and gluon 
condensate is important in the dynamical hQCD model to produce the correct light flavor 
meson spectra. 

The pion form factor and the vector couplings are also investigated in the dynamical 
hQCD model. In this part, the results in our model is not worse than the other models. 

However, we emphasize that we have offered a systematic framework to describe the 
gluodynamics and chiral dynamics. The input in our model is basically the gluon conden- 
sate represented by the chiral condensate a, and a current quark mass ruq, which are 
the same as in the soft-wall model. Just solve the deformed warp factor self-consistently. 
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we can produce the glueball spectra, the hnear heavy quark potential as well as light flavor 
meson spectra in very well agreement with lattice and experimental data. 
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